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We study the qualitative properties of the class of spatially homogeneous Bianchi VIo cosmological models 
containing a perfect fluid with a linear equation of state, a scalar field with an exponential potential and a 
uniform cosmic magnetic field, using dynamical systems techniques. We find that all models evolve away from 
an expanding massless scalar field model in which the matter and the magnetic field are negligible dynamically. 
We also find that for a particular range of parameter values the models evolve towards the usual power-law 
inflationary model (with no magnetic field) and, furthermore, we conclude that inflation is not fundamentally 
affected by the presence of a uniform primordial magnetic field. We investigate the physical properties of the 
Bianchi I magnetic field models in some detail. 

I. INTRODUCTION 

There are many observations that imply the existence of magnetic fields in the Universe, and recent observations have produced 
firmer estimates of the strength of magnetic fields in interstellar and intergalactic space (Kronberg, 1994; Kronberg et al., 1992; 
Wolf et al., 1992). The observations to date place an upper bound on the strength of a cosmic magnetic field, but are not 
conclusive as to whether such a field exists (see Vallee 1990a,b). Although cosmologists have investigated the possible existence 
of a homogeneous intergalactic magnetic field of primordial origin both from a theoretical and an observational point of view for 
many years, research on cosmological magnetic fields has been rather marginal despite their potential importance. The reasons 
could be the perceived weakness of the field effects or the lack, as yet, of a consistent theory explaining the origin of cosmic 
magnetism. However, this situation has changed considerably recently (cf. Grasso and Rubenstein, 2000). 

Primordial magnetic fields introduce new ingredients into the standard picture of the early Universe. Such a field would affect 
the temperature distribution of the microwave background radiation, primeval nucleosynthesis and galaxy formation. But its 
most direct observational effect is the Faraday rotation it would cause in linearly polarized radiation from observed extragalatic 
radio sources. Fundamental properties of magnetic fields include their vectorial nature, which inevitably couples the field to the 
spacetime geometry, and the resulting tension (i.e., the negative pressure) exerted along the field's lines of force (Parker, 1979; 
Matravers and Tsagas, 2000 ). The implications of such an interaction are both kinematical and dynamical; kinematically, the 
magneto-curvature effect tends to accelerate positively curved perturbed regions, while it decelerates regions with negative local 
curvature and, dynamically, the most important magneto-curvature effect is that it can reverse the pure magnetic effect on density 
perturbations (Tsagas and Barrow, 1997 and 1998; Tsagas and Maartens, 2000 and 2001). 

In a recent analysis (Matravers and Tsagas, 2000) the kinematics were considered and, assuming a spacetime filled with 
a perfectly conducting barotropic fluid permeated by a weak primordial magnetic field while treating the energy density and 
the anisotropic pressure of the field as first-order perturbations upon the Friedmann-Lemaitre-Robertson-Walker (FLRW) back- 
ground, it was shown how a cosmological magnetic field can modify the expansion rate of an almost-FLRW universe. The effects 
of the interaction between magnetism and geometry in cosmology are subtle, and it was argued that the magneto-curvature cou- 
pling could make the field into a key player irrespective of the magnetic strength. Indeed, even weak magnetic fields can lead 
to appreciable effects, provided that there is a strong curvature contribution. This was illustrated by studying spatially open cos- 
mological models containing 'matter' with negative pressure, and it was found that the phase of accelerated expansion, which 
otherwise would have been inevitable, may not even happen. This leads to the question of the efficiency of inflationary models 
in the presence of primordial magnetism. Recall that an initial curvature era was never considered as a problem for inflation, 
given the smoothing power of the accelerated expansion. However, this may not be the case when a magnetic field is present, 
no matter how weak the latter is. In more recent work Tsagas (2000) has further studied the vector nature of magnetic fields and 
discussed the unique coupling between magnetism and spacetime curvature in general relativity which gives rise to a variety of 
effects with important implications. 

Observations of the high degree of isotropy of the cosmic microwave background indicate that the Universe is 'almost' 
isotropic and spatially homogeneous (at least since the time of last scattering). Theoretical support for this belief comes from the 
so-called Ehlers-Geren-Sachs (1964) Theorem. Clarkson and Coley (2001) have proven a generalization of this theorem and have 
consequently shown that any strong magnetic fields in the Universe are ruled out. This theoretical result is model-independent 
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and includes the case of inhomogeneous magnetic fields. In further work (Clarkson et al., 2001) numerical constraints are placed 
on all types of primordial and protogalactic magnetic fields in the Universe from cosmic microwave background data. 

We shall study open spatially homogeneous cosmological models containing both a uniform magnetic field and a scalar field 
here, partially in an attempt to address some of the questions raised above. Scalar field cosmological models are of great 
importance in the study of the early Universe. Models with a variety of self-interaction potentials have been studied, and one 
potential that is commonly investigated and which arises in a number of physical situations has an exponential dependence on 
the scalar field (Wetterich, 1988; Ferreira and Joyce, 1997). There have been a number of studies of spatially homogeneous 
scalar field cosmological models with an exponential potential, with particular emphasis on the possible existence of inflation in 
such models (Halliwell, 1987; Kitada and Maeda, 1993; and Billyard et al., 1999, and references within). 

Dynamical systems methods for analysing the qualitative properties of cosmological models have proven very useful (Wain- 
wright and Ellis, 1997; Coley, 1999). A qualitative analysis of cosmological models with matter and a uniform magnetic field 
has been presented previously (Collins, 1972; LeBlanc et al., 1995). A universe with a primordial magnetic field is necessarily 
anisotropic. Thus, in order to investigate the influence of the magnetic field on the dynamics of the universe the Einstein field 
equations must be analysed in models more general than the FLRW models. The simple classes of the Bianchi I and Kantowski- 
Sachs models were discussed (cf. Jacobs, 1969), and Collins (1972) applied techniques from the theory of planar dynamical 
systems to prove qualitative results about the evolution of the class of axisymmetric Bianchi I cosmologies with matter and a 
primordial magnetic field. 

In LeBlanc et al. (1995) the Einstein-Maxwell field equations for orthogonal Bianchi VIo cosmologies with a 7-law perfect 
fluid and a pure, homogeneous source-free magnetic field were written as an autonomous differential equation in terms of 
expansion-normalized variables. It was shown that the physical region of state space is compact, and that the differential 
equation admits certain invariant sets and monotone functions which play an important role in the analysis. A complete analysis 
of the stability properties of the equilibrium points of the differential equation and a description of the bifurcations that occur as 
the equation of state parameter 7 varies was given . The associated dynamical system was studied and the past, intermediate and 
future evolution of these models was determined. All asymptotic states of the models, and the likelihood that they will occur, 
were described. In particular, oscillatory behaviour also occurs in cosmological models with a magnetic field (and in Einstein- 
Yang-Mills theory in general). Further work on spatially homogeneous models with a magnetic field and a non-tilted perfect 
fluid has been carried out recently (LeBlanc, 1997 and 1998). In particular, Weaver (2000) has generalized (to the non-polarized 
solutions) the work of Leblanc et al. (1995), and rigorously shown that the evolution toward the singularity is oscillatory in 
Bianchi VI vacuum models. 

The reasons for concentrating on the Bianchi VI models are ones primarily of mathematical simplicity. It is well known 
(Hughston and Jacobs 1970) that a pure magnetic field is only possible in Bianchi cosmologies of types I, II, VIo, VIIo (in class 
A) and type III (in class B). For types VI and VII , the algebraic constraints that arise from the Einstein field equations imply 
that the shear eigenframe is Fermi-propagated, which in turn implies that the remaining field equations reduce to an autonomous 
differential equation with a polynomial vector field. Finally, for type VIo but not for VIIo, the physical region of state space is 
compact. We note that since a magnetic field is not compatible with Bianchi types VIII and IX (for example), the class of Bianchi 
VIo magnetic cosmologies is of the same generality as the Bianchi VIII/IX models (without magnetic field) (Wainwright 2000). 

In this paper we shall investigate the class of Bianchi VI models with barotropic matter, a scalar field with an exponential 
potential and a uniform magnetic field. We shall discuss some general properties of the complete class of models, and then 
investigate the special case of the Bianchi I system in detail. We note that this is not the general class of Bianchi I cosmological 
models; we are only considering those Bianchi I models that occur as subcases of the Bianchi VI models with one magnetic 
degree of freedom. It is also of interest to study the Bianchi II cosmologies, since although they are very special within the full 
Bianchi class, they play a central role since the Bianchi II state space is part of the boundary of the state space for all higher 
Bianchi types (i.e., all types except for I and V). The Bianchi II models will be investigated in detail elsewhere (Quinlan, 2002). 

The Bianchi I system in the absence of a scalar field was first qualitatively analysed in (Collins, 1972). A qualitative analysis 
of the Bianchi I models in the absence of a magnetic field was given in (Billyard et al., 1999); in this work the well-known 
power-law inflationary solution (Wetterich, 1988; Kitada and Maeda, 1993) was shown to be a stable attractor for an appropriate 
parameter range in the presence of a barotropic fluid in all Bianchi class B models (previous analysis had shown that this power- 
law inflationary solution is a global attractor in spatially homogeneous models in the absence of a perfect fluid, except for a 
subclass of Bianchi type IX models which recollapse). One of the aims of the present analysis is to study the stability of this 
model with respect to magnetic field perturbations. 



II. THE BIANCHI VIq MODELS 

We shall follow the approach of LeBlanc et al. (1995) in which the theory of dynamical systems was used to give a detailed 
analysis of the evolution of orthogonal Bianchi cosmologies of type VIo with a perfect fluid and a magnetic field as source. This 
work extended that of Wainwright and Hsu (1989) who had studied perfect fluid Bianchi cosmologies of class A using state 
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variables that are dimensionless and have a direct physical or geometric interpretation in terms of the shear of the cosmological 
fluid, the spatial curvature and the magnetic field, leading to a state space that is a compact subset of M 5 , which implies that 
a unified treatment of the asymptotic behaviour of the models at early and late times can be given. We note that all of the 
equilibrium points of the differential equation correspond to self-similar exact solutions of the Einstein field equations. 

An invariant orthonormal frame of vector fields on the spacetime are introduced in which one vector is aligned along the 
fluid flow vector so that the remaining three spatial vectors (triad) span the tangent space orthogonal to the fluid flow at each 
point of the group orbits. The commutation functions of this frame are then taken as the basic variables. The Einstein-Maxwell 
field equations for a pure magnetic field are then derived in the orthonormal frame formalism, giving rise to a set of evolution 
equations for the shear variables (<7 a p), the curvature variables (n a/ g), the energy density (conservation equation), the magnetic 
field (Maxwell equation) and a first integral (the Friedmann constraint), where the remaining freedom in the choice of spatial 
tetrad was used to simultaneously diagonalize a a p and n a p. 

Expansion-normalized variables are then introduced, leading to a reduced set of evolution equations (the dynamical system) 
for the dimensionless shear variables E±, the spatial curvature variables N a , the dimensionless magnetic field variables and the 
density parameter (defined as usual). Defining a dimensionless time variable, the reduced evolution equations in the special case 
of a Bianchi VIo were then established (in this case there is a single dimensionless magnetic field variable, M). 



A. The Equations 

We assume a scalar field with an exponential potential V = Voe , where Vb and k are positive constants, and a separately 
conserved perfect which satisfies the barotropic equation of state p = (7 — l)p, where the constant 7 satisfies < 7 < 2 
(although we shall only be interested in the range < 7 < 2 here). The energy conservation and Klein-Gordon equations 
become: 

p = -3jHp, (2.1) 
= -3H<f> - kV, (2.2) 

where H is the Hubble parameter, an overdot denotes ordinary differentiation with respect to time t, and units have been chosen 
sothat87rG = 1. 

We also assume a uniform magnetic field (in the 'x' direction) which satisfies the Maxwell equations. The form of this 
magnetic field, which we shall denote here as M 2 , is given in Collins (1972) or LeBlanc et al. (1995), but its precise form is 
unimportant since we use the generalized Friedmann constraint to eliminate M 2 from the evolution equations. 

We define 

Q = - $ = - * = — (2 3) 

3ff 2 ' 3ff 2 ' ( ' 

two normalized (dimensionless) shear variables X+ and £_ (appropriately combinations of the derivatives of the metric func- 
tions divided by the Hubble parameter; Wainwright and Ellis, 1997)), and the new logarithmic time variable r by 

d i=H. (2.4) 

at 

The evolution equations for the quantities 

X = (£+,£_,$,4<,0,A + ,A_) e K 7 , (2.5) 



are then as follows: 



E' + = (q - 2)S+ + 2(1 - A^ 2 - $ 2 - * - Y? + - Y?_ - 0) - 2A 2 , (2.6) 

=(q- 2)S_ - 2A+A_, (2.7) 

$' = (4-2)$-^*, (2.8) 

= (2^ + 2 + V6fc$)*, (2.9) 

n' = (2<7-3 7 + 2)ft, (2.10) 

N' + = (g + 2£ + )A + +6£_A_, (2.11) 

N'_ = (q + 2£+)A_ + 2£_A + . (2.12) 
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where a prime denotes differentiation with respect to the time r and the deceleration parameter q is defined by q = —(l+H'/H), 
where 



q = 1 - N 2 _ + Y? + + + i(3 7 - 4)fi + $ 2 - 2*. (2.13) 

The physical state space is restricted to il > and N 2 < 3N 2 (where we include equality to include models of Bianchi types I 
and II; LeBlanc et al., 1995). Since N- = implies N'_ = 0, the subsets N- > and N- < are invariant. Due to a discrete 
symmetry in evolution equations we can, without loss of generality, restrict attention to the subset iV_ < 0. 
The magnetic field, defined through the first integral, is given by 

|iW 2 = 1 -N 2 _ - $ 2 - Y? + - Y, 2 _ - n, (2.14) 

and since M 2 > 0, all of the variables are bounded: 

< {Z 2 + ,Z 2 _,<i> 2 ,-$,n,Nl/3,-N-} < 1. (2.15) 

Without loss of generality, we can restrict attention to the subset M > 0. There is also an auxiliary equation for the magnetic 
field: 

M 1 = (q - 1 - 2S+)M. (2.16) 



B. Discussion 



A number of invariant sets of the physical state space can be identified (cf., LeBlanc et al., 1995, and Billyard et al., 1999). In 
addition, there are a number of monotonic functions that exist in various invariant sets (particularly on the boundary; cf. p. 521 
LeBlanc et al., 1995). However, we shall not present them here. Indeed, there are many (over fifty) equilibrium points, the vast 
majority of which are saddles, and hence in this section we shall simply comment on some of the more physically interesting 
local properties of the models. In the next section we shall investigate the Bianchi I models more rigorously, and a monotonic 
function in this invariant set will be given explicitly, in order to discuss some aspects of the intermediate or transient behaviour 
of the models. 

• The equilibrium point V : (£+, £_, Q, N + , iV_) = (0,0, 1- ^,0, 0,0), M = corresponds to the power-law 
inflationary zero-curvature (flat) FLRW model. The eigenvalues of V are {fc 2 -4, A: 2 -37, ±(/c 2 -2)(x2), ±(fc 2 -6)(x3)}. 
V exists for k 2 < 6 and is a sink for k 2 < 4 and k 2 < 3j, and is inflationary (q < 0) for k 2 < 2. In fact it is the global 
attractor for this range of parameter values; this result generalizes previous work by including a cosmic magnetic field (cf. 
Billyard et al., 1999). 

• We are particularly interested in whether there are any sinks (or sources) with M ^ in the Bianchi VIo state space. 
We know from LeBlanc et al. (1995) that there are no attractors in the absence of a scalar field, so hence we must have 
M 0, <3> 2 ^ (and ^ 0, since this leads to M = 0). From equation (2.10) we see that O approaches zero at late 
times, so we also take = (and q < 2). From equations (2.8), (2.9) and (2.16) we have that (for M ^ 0) 

$=^(1 + S+), * = _L ( i + E+)( i_ 2 £ + ); g =l + 2S + . (2.17) 

We are particularly interested in the inflationary case in which q < 0; i.e., —1 < E + < — |. If N-,N + (S-) are 
not zero, equations (2.11) and (2.12) yield N + = ±V37V_, £_ = T^g(l + 4S+), whence equation (2.7) then yields 

N 2 = (1 — 2£ + )(l + 4S + ), which is negative for the range of E + under consideration! (There do exist equilibrium 
points with N- ^ 0, N + ^ and £_ ^ 0, but these all have q > (i.e., are non-inflationary) and are saddles.) Hence we 
consider N_ = N + = £_ = 0. Equations (2.6) and (2.13) then yield (since £+ ^ \), £+ = -1 + 3fc 2 /(8 + k 2 ). This 
equilibrium point, denoted V in section V, is physical and in the state space when 4 < k 2 < 8 and is always a saddle (N+ 
is increasing as orbits evolve away from V to the future). However, since q < only for k 2 < | (S + < — 5), this point is 
never inflationary. 

We conclude that inflation is not fundamentally affected by the presence of primordial magnetism in the models under study 
here. 
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• In LeBlanc et al. (1995) it was shown that in the absence of a scalar field that for fi>0, 1<7<2 the global sink 
is the equilibrium point PM(VIo) so that at late times all such magnetic Bianchi VIo cosmologies are approximated by 
the dynamics of a self-similar Dunn-Tupper (1980) magneto-vacuum Bianchi VIo model. This equilibrium point has an 
analogue here with $ = ^ = 0, but is now a saddle. Indeed, it can be easily shown that none of the equilibrium points 
with vanishing scalar field can be sinks. In particular, the sink in the case 7 = 1 in LeBlanc et al. (1995) (subset of 
£M(VI )) becomes a saddle. 

• All equilibrium points corresponding to models with a non-zero matter fluid, including those corresponding to the flat 
FLRW perfect fluid solution (Einstein de Sitter) and the flat FLRW scaling solution, are saddles. 

In the above we have made some passing comments concerning the intermediate behaviour of the models, but we have not 
discussed this in detail. In the next section we will discuss the intermediate behaviour of the (subset) of Bianchi I models. The 
transient behaviour is important for assessing the physical significance of the models. For example, LeBlanc et al. (1995) also 
considered whether any of the Bianchi VIo magnetic cosmologies are compatible with observations of a highly isotropic uni- 
verse. For a subset of the Bianchi VI cosmologies, the orbits will initially be attracted to the equilibrium point T corresponding 
to the flat FLRW model, but will eventually be repelled by it. For such a model there will be an epoch during which the model is 
close to isotropy and hence will be compatible with observations. This is the phenomenon of intermediate isotropization. It will 
not occur for a randomly selected magnetic Bianchi VI model, but there will be a finite probability that an arbitrarily selected 
model will undergo intermediate isotropization (the probability will depend on the desired closeness to isotropy and the length 
of the epoch of isotropization). 

• There is an equilibrium set JC: (£+, £_, 0, 0, 0, 0) where 1 - T,\ - IP - $ 2 = (and M = 0, q = 2). (These points 
reduce to the equilibrium set JCm in the absence of a magnetic field consisting of Jacobs type I solutions with no matter and 
a massless scalar field - Billyard et al., 1995). The eigenvalues are {0, 0, 2(1 - 2S+), 3(2 - 7), 2(1 + S+± V3X_), 6(1 + 
^|)}- There are two zero eigenvalues (corresponding to the fact that JC is a 2-parameter family of equilibrium points). 
An analysis shows that there is always a subset of JC that acts as sources (these are, in fact, global sources); all others are 
saddles. 

This contrasts with the situation in the absence of a scalar field in which LeBlanc et al. (1995) showed that there are no 
equilibrium points that are sources and argue that generically an orbit with fi > will pass through a transient stage and then 
approach the Kasner circle into the past (consisting of Kasner, Bianchi type I vacuum models). Since these equilibrium points 
are saddles, the orbit subsequently leaves along a uniquely determined heteroclinic orbit and then approaches the Kasner circle 
again. This process then continues indefinitely, and the orbit follows an infinite heteroclinic sequence of Rosen orbits and Taub 
orbits joining Kasner equilibrium points undergoing Mixmaster-like chaotic oscillations (the past attractor here consists of the 
union of the Kasner circle and a family of Rosen magneto-vacuum type I orbits and Taub vacuum type II orbits - see also 
Weaver, 2000). 



III. BIANCHI I MODELS 



The evolution equations for the Bianchi I models are obtained by restricting the above equations to the zero-curvature case 
(i.e., N± — 0), whence we obtain 

X = (£+,£_, $, <J/, O) e M 5 , (3.1) 

and the evolution equations are as follows: 

* -T, 2 + -Y, 2 _-Q), (3.2) 

(3.3) 

(3.4) 

(3.5) 
(3.6) 

where a prime denotes differentiation with respect to the time r and the deceleration parameter q is now given by 

q = 1 + Z 2 + + + i(3 7 -4)fi + $ 2 - 2*. (3.7) 





= (?- 


2)S 4 


. + 2(l-$ : 
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= (?- 
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= (2q- 
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= (2«- 


-3 7 - 


f2)Q, 
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The magnetic field satisfies 



^M 2 = 1 - $ 2 - * - E 2 . - X 2 - fl, (3.8) 

and since M 2 > 0, all of the variables are bounded: < {X^, X 2 , $ 2 , ft} < 1. 

The equations (3.1)-(3.6) can also be easily derived from the metric directly in coordinate form - see the Appendix. 



A. Invariant sets and Monotonic functions 



The existence of monotonic functions rule out periodic and recurrent orbits in the phase space and enables global results to be 
obtained. When O ^ 0, we can define Z = ftXl 2 , whence from the evolution equations above we find that 

Z' = 3(2 - j)Z, (3.9) 

so that Z is a monotonically increasing function for 7 < 2 (LeBlanc et al., 1995). From the monotonicity principle (Wainwright 
and Ellis, 1997) we can deduce that fl — > at early times and X 2 . — > at late times. 
When ft — 0, we note that 

< ? -2 = 3(M 2 + X 2 h + X 2 _+$ 2 -l) <0, (3.10) 

and so from equation (3.3) the variable X_ is itself a monotonic function. 

There are a number of invariant sets. M = is an invariant set corresponding to the absence of a magnetic field. The dynamics 
in this invariant set was studied in Billyard et al. (1999). $ = * = is an invariant set corresponding to the absence of a scalar 
field. The dynamics in this invariant set was studied in Collins (1972). ft = is an invariant set corresponding to the absence 
of barotropic matter. This invariant set is important for studying the the early time asymptotic behaviour. Finally, X_ = is an 
invariant set which is important in the study of late time asymptotic behaviour. 



B. Stability of Equilibria 

Let us present all of the equilibrium points and their associated eigenvalues. We shall assume here that 1 < 7 < 2. The case of 
the bifurcation value 7 = 2 will be discussed in Quinlan (2002). All equilibrium points correspond to self-similar cosmological 
models. 

Equilibrium Points: {X + ,X_,$,*,f2}: 



1. Zero magnetic field 



V - Power-law inflation (flat FLRW): 



k_ fc 2 
V6' "6 



0,0, — ^=,1-^,0J>, M = 0, q=\k 2 -l. 



Eigenvalues: [fc 2 - 4, fc 2 - 3 7 , ifc 2 - 3, ifc 2 - 3, ifc 2 - 3] 

Exists provided fc 2 < 6. Sink if fc 2 < 4 and fc 2 < 37; saddle otherwise. Inflationary if fc 2 < 2 (note that this is the only 
equilibrium point that can be inflationary). 

T s - Matter-scaling (flat FLRW; Copeland et al., 1998): 



Eigenvalues: [|( 7 - 2), |( 7 - 2), 3 T - 4, § (7 - 2) ± ^[2 - -,) (24 7 ^ + P [2 - 97]) 
Exists provided fc 2 > 37. Sink if 7 < |; saddle otherwise. 
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• T - Einstein-de-Sitter (flat FLRW): 

{0, 0, 0,0, 1} M = 0, q = §7 - 1 > 0. 



Eigenvalues: [§ ( 7 - 2) , § ( 7 - 2) , § ( 7 - 2) , 3 7 , 3 7 - 4] 
Saddle. 

• JC - Kasner vacuum/massless scalar field: 



{3 



$ 2 = 1, * = n = 0} , M = 0, g = 2. 



Eigenvalues: [0, 0, 3(2 - 7 ), 2 - 4E+, 6 + V6fc$] 

This is a 2 -parameter family of equilibrium points. Sources when E + < i; saddles otherwise. 

2. Non-zero magnetic field 

These all correspond to anisotropic models + E 2 . ^ 0). 
• J - Matter and no scalar field (Jacobs magnetic field model; Jacobs, 1969): 

{i(3 7 -4),0,0,0,|(l-i 7 )}, Af 2 = ±(2- 7 )(3 7 -4), q = | 7 - 1 > 0. 



( 7 - 2) , | ( 7 - 2) , 3 7 , | {( 7 - 2) ± V(2- 7 )(3 7 2 - 17 7 + 18)} 



Eigenvalues: 

Exists provided 7 > | (same as J 7 for 7 = |). Saddle. 
• 5 - Matter and scalar field: 

|i(3 7 -4),0,-^J,^(2- 7 ),^(4fc 2 -fc 2 7 -8 7 )|, Af 2 = ±(2- 7 )(3 7 -4), ? =§ 7 -l>0. 



Eigenvalues: 



( 7 -2),|{2( 7 -2)±I v /2(2- 7 ) (a±b)} 



where a = 24 7 2 + 3 7 2 fc 2 - 26fc 2 7 + 20fc 2 , and b = v/(18fc 2 7 - 4fc 2 + a - 32fc 7 ) (18fc 2 7 - 4fc 2 + a + 32fc 7 ). 

4fc 2 



Exists provided 7 > | (same as for 7 = |) and 7 < p^g . Sink. 



• V - Scalar field and no matter: 



fc 2 -4 



,0,- 



2\/6fc 



12- 



8 — fc 2 



fc 2 + 8' ' fc 2 + 8' (p + 



r,0^, AT 2 = 2 



(fc 2 - 4) (8 - fc 2 ) 
(P + 8) 



5fc 2 



fc 2 + ; 



> 0. 



Eigenvalues: 



o «~ —a q fc — 8 o* 



fc 2 +8 



± 



17fc 2 -72 
fc 2 -8 



Exists provided 4 < fc 2 < 8. Sink when 7 > i fc g ; saddle otherwise 



fc 2 +8' 



C. Discussion 

In Table | we summarize the equilibrium points and their stability for different values of the parameters ( 7 , fc 2 ). In Figure [l] 
we show the sinks in (7, fc 2 )-parameter space. The equilibrium point V represents the power-law FLRW solution with no matter 
and no magnetic field, and is inflationary for fc 2 < 2. It is a global attractor for fc 2 < 3 7 in the Bianchi I invariant set. Note that 
V is a sink in the full class of Bianchi VIq models, unlike all other sinks in Figure 1 which are only attractors in the Bianchi I 
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invariant set. Also note that the monotonic function Z in the Bianchi I invariant set indicates that O — > at early times (at the 
sources) and £_ — > at late times (at the sinks), consistent with the results in Table |. 

The equilibrium points S and V, in which all of the matter fields are non vanishing and in which matter field (only) vanishes, 
respectively, correspond to new exact self-similar magnetic field cosmological models. The new solutions are given explicitly in 
the Appendix. 

Heteroclinic sequences (and hence the transient behaviour of the Bianchi I models) can easily be deduced from the Table. As 
a particular example, for 7 = 5/3 and 5 < k 2 < 40/7, a subset of K, constitute sources (the remainder are saddles), V is a 
sink, J is a saddle (that lies in the $ = 0, $ = (M = 0) submanifold) and T, V and Ts are saddles (that lie in the isotropic 
£+ = 0, £_ = (M = 0) submanifold-in this submanifold JC is a source, T and V are saddles, and Ts is a sink). The 
heteroclinic sequence in this case is presented in Figure |[ 

From the the Table and Figures we can deduce both the asymptotic and the intermediate behaviour of the Bianchi I models. 
Note that as 7 increases, the magnetic field becomes increasingly important at late times, consistent with the observation of 
Wainwright (2000) (in the absence of a scalar field). 



(a) 1 < 7 < s 



k 2 


37 4 6 8 


V 


SINK 


SADDLE 


DOES NOT EXIST 


T s 


DOES NOT EXIST 


SINK 


T 


SADDLE 


V 


DOES NOT EXIST 


SADDLE 


DOES NOT EXIST 



(b) § < 7 < 2 



k 2 


4 3 7 6 87/(4-7) 8 


V 


SINK 


SADDLE 


DOES NOT EXIST 


Ts 


DOES NOT EXIST 


SADDLE 


T 


SADDLE 


J 


SADDLE 


S 


DOES NOT EXIST 


SINK 


V 


DOES NOT EXIST 


SINK 


SADDLE 


DOES NOT EXIST 



TABLE I. Summary of all of the equilibrium points (as defined in the text) and their nature for different {7, fc 2 }-parameter values in the 
Bianchi I models with (a) 1 < 7 < | and (b) 7 > |. Note that for | < 7 < 2, 4 < 37 < 87(4 — 7)" 1 < 8, and in (b) we have identified 
87(4 — 7) -1 > 6 (i.e., 7 > ^) to be specific. For all values of 7, a subset of JC are sources (the remainder are saddles). 



IV. CONCLUSIONS 

We have investigated the class of open Bianchi VIo universe models with barotropic matter, a scalar field with an exponential 
potential and a uniform magnetic field. We discussed some of the general properties of the models by utilyzing dynamical 
systems techniques. The equilibrium point V (with M = 0), corresponding to the power-law inflationary flat FLRW model, is a 
global sink for k 2 < 2. Hence all models are future aymptotic to this inflationary attractor for these parameter values. There is 
an equilibrium set JC with 1 — Y?L — — $ 2 = (and M = 0, q = 2), a subset of which are global sources. Hence all models 
are past asymptotic to massless scalar field models with no matter and no magnetic field. 

A (partial) analysis of the saddles was undertaken in order to determine some of the transient features of the models. In 
particular, we found that there are no equilibrium points with M 7^ in the Bianchi VIo state space which are inflationary, 
and hence concluded that inflation is not fundamentally affected by the presence of a uniform primordial magnetic field in these 
models. 

This latter result is not necessarily inconsistent with the conclusions of Matravers and Tsagas (2000), since only a uniform 
magnetic field was considered here. Indeed, it is the non-uniform magnetic field gradients that give rise to the magneto-curvature 
effects in their work which modify the cosmological expansion rate of an almost-FLRW universe and that can have undesirable 
implications for inflationary models (and may even prevent inflation taking place in the presence of primordial magnetism). It 
might be thought that on large scales that the magnetic field will be approximately homogeneous, but Matravers and Tsagas 
(2000) have argued that even weak magetic effects may have significant cosmological consequences. The drawback of their 
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work is that only a local perturbation analysis was performed and questions of genericity and long term behaviour cannot be 
easily addressed. Clearly, an investigation of the qualitative properties of a class of scalar field cosmological models with an 
inhomogeneous magnetic field would extend and generalize both the work of Matravers and Tsagas (2000) and the (restricted) 
analysis of a uniform magnetic field here. 

In order to investigate the intermediate behaviour of the models, and hence their physical properties, we discussed the (subset) 
of Bianchi I models in more detail. We found the equilibrium points S and V, in which all of the matter fields are non vanishing 
and in which the matter field (only) vanishes, respectively (corresponding to new exact self-similar magnetic field cosmological 
models), act as sinks in the Bianchi I subset. Heteroclinic sequences (and hence the transient behaviour of the Bianchi I models) 
were discussed. We found that as 7 increases, the magnetic field becomes increasingly important at late times. 
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APPENDIX A: 

In the main text we followed the dynamical systems theory approach of LeBlanc et al. (1995) to analyse the evolution of 
orthogonal Bianchi cosmologies of type VIo with a uniform magnetic field. This approach uses an invariant orthonormal frame 
in which the commutation functions are the basic variables. Expansion-normalized variables are then introduced, and a reduced 
set of evolution equations (the dynamical system) for the dimensionless shear variables and spatial curvature variables as well 
as the dimensionless magnetic field variables and the density parameter are then obtained. All of the equations in the text follow 
from the analysis of LeBlanc et al. (1995), wherein all quantities are properly defined. However, a coordinate approach is also 
possible, and it may be of use for the reader to get a better sense of how the equations are derived, at least in the case of the 
Bianchi I models. 

The metric for the class of anisotropic Bianchi I models, which are the simplest spatially homogeneous generalizations of the 
flat FLRW models which have non-zero shear but zero three-curvature, is given by 

ds 2 = -dt 2 + X 2 dx 2 + Y 2 dy 2 + Z 2 dz 2 , (Al) 

where X, Y, Z are functions of t only. The expansion rate is given by 

^3tf=§4 + §, (A2) 
and when Y = Z (the LRS subcase) there is only one independent rate of shear, which is given by 

(A3) 




Normalized (dimensionless) variables and a logarithmic time variable r are defined by equations (2.3) and (2.4), which then lead 
to the evolution equations (3.1)-(3.6) in coordinate form (Collins, 1972). 

1. Exact Solutions 

We can write down the exact solutions corresponding to the equilibrium points S and V in the text. 

- S: Since q = §7 — 1, we can immediately integrate the Raychaudhuri equation H' / H = —(1 + q) to obtain H = 3^t _1 . 
Since £_ = and £+ = \(3>y - 4), it follows that Y = Z, so that a = S + x H = \{l - j-)^ 1 . The metric is 
consequently of power-law form, viz., 



ds 2 = -dt 2 + t 2pi dx 2 + t 2p2 (dy 2 + dz 2 ), (A4) 
where pi and p 2 are constants. From the expressions for H and a we then find that pi + 2p 2 = ~ and p 2 — pi — ' 



and hence 



7 """^ ft — 2 7 
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2 1 

Pi = - - 1, P2 = ■ (A5) 

7 2 

The geometry is that of the Jacobs magnetic field model (Jacobs, 1969), and the only difference here is that the matter 
fields constituting the source are given by 



P = XlH>=( 4 #-*>-*')t->, Hl ^MH=^-^- A \-\ (A6) 



and the scalar field is given by 



where d>n = 



fa = In ^ ^fc^vb" 1 ) ■ This solution exists provided 7 > 



- V: In this solution q = 5 fc fc a~ +8 8 , so that = 1 . £_ = and E + = 2pr|, and so the metric is of the form of 

equation (A4), where now 



P! = ^, P 2 = i <A8, 



There is no perfect fluid in this solution (p = 0), and the magnetic field is given by 



Hi = jHEMl r i. ( A9) 

The scalar field is given by equation (A7) with O = \ In ( fr^r ) ■ This new solution exists provided 4 < k 2 < 8. 
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